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Axion-like particles (ALPs) and photons inter-convert in the presence of a magnetic field. At keV 
energies in the environment of galaxy clusters, the conversion probability can become unsuppressed 
for light ALPs. Conversion of thermal X-ray photons into ALPs can introduce a step-like feature 
into the cluster thermal bremsstrahlung spectrum, and we argue that existing X-ray data on galaxy 
clusters should be sufficient to extend bounds on ALPs in the low-mass region ma 1 x 10“^^ eV 
down to M ~ 7 X 10^^ GeV, and that for 10^^ GeV < M < 10^^ GeV light ALPs give rise to 
interesting and unique observational signatures that may be probed by existing and upcoming X- 
ray (and potentially X-ray polarisation) observations of galaxy clusters. 


Galaxy clusters are the largest gravitationally bound 
objects in the universe and provide a powerful testing 
ground for theories of new physics. Recently, it has been 
appreciated that clusters are highly efficient at inter¬ 
converting light axion-like particles (ALPs) and photons 
[1-3] (see [4-6] for some earlier work). In this paper, we 
show that the absence of large distortions of the cluster 
thermal X-ray bremsstrahlung spectrum may be used to 
derive the strongest bounds to date on the ALP-photon 
coupling for light ALPs. 

ALPs, reviewed in [7], arise in many theories of physics 
beyond the Standard Model and are ubiquitous in string 
theory compactifications. The phenomenological low- 
energy Lagrangian for ALPs and photons is given by, 

£ = Id.ad^a + , ( 1 ) 


where M is the ALP-photon coupling, and rria denotes 
the mass of the ALP. 

Significant observational and theoretical effort has 
gone into searching for and constraining such parti¬ 
cles. The strongest current bound arises from obser¬ 
vations of SN1987A, leading to M > 2 x 10“ GeV [8- 
10]. Planned experiments such as ALPS-II or lAXO 
are expected to produce competitive bounds [11, 12], 
e.g. M > 3 X 10^^ GeV in the lAXO best-case-scenario. 

For massive ALPs, a recent paper used the absence of 
CMB distortions through clusters to produce bounds on 
M [13], although as this relies on resonance effects it is 
only relevant for a small range of ALP masses. 

Here, we argue that current and future data from X-ray 
observations of galaxy clusters can significantly improve 
the bounds on M for the entire small nia region (< 1 x 
10“^^ eV), reaching down to around M ~ 10^^ GeV. 


Photon-ALP conversion in galaxy clusters 


As is well known, the third term of the Lagrangian (1) 
induces ALP-photon inter-conversion in the presence of 
coherent magnetic fields [14-16]. The linearised equa¬ 
tions of motion for the modes of energy w propagating in 
the z-direction are given by. 
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for the ALP-state |a) and the two transverse polarisations 
of the photon beam, [ya,) and \Uy)- Here Aj = 

Ajai = Bi/2M and A^ = —m„/2a;. The plasma fre¬ 
quency of the surrounding medium is given by, 
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( 3 ) 

Ap accounts for Faraday rotation. However this effect 
depends on wavelength as and thus is negligible at 
X-ray energies. 

For a constant magnetic field in a domain of length 
L, the photon-to-ALP conversion probability for unpo¬ 
larised light is given by, 
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Here B± denotes the component of the magnetic field 
that is perpendicular to the ALP wave vector. The factor 
of I accounts for the fact that only one polarisation state 
of light participates in the mixing. 

It is instructive to consider the typical values of 0 and 
A in galaxy clusters. For ma = 0 we have, 
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For 0, A <C 1, conversion is quadratic in both size and 
coherence length of the magnetic field, 
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and ALP-photon conversion is energy-independent. 
However, for 0 <C 1 and A ^ 1, the conversion prob¬ 
ability is energy-dependent and progressively suppressed 
at lower energies. 
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We can define a critical energy marking the crossover 
between the regimes, 
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and thus for typical clusters and for M > lO^^GeV, the 
crossover between regimes occurs within the range of ob¬ 
served X-ray photon energies. We now show that this can 
lead to substantial distortions of the thermal spectrum of 
the intra-cluster medium, using the Coma cluster as our 
prime example. 
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X-ray emission from galaxy clusters 


matter line reported in [20] is observed as a per cent level 
effect above the background, thereby requiring the abil¬ 
ity to characterise the background thermal emission to 
the same level of accuracy. 

Thus, galaxy clusters provide an intrinsic, diffuse, 
bright and well characterised source of X-ray photons. 


Cluster spectral distortions from ALPs 

If a significant proportion of the thermal photons were 
to convert to ALPs, the resulting X-ray spectrum would 
be distorted. The photon-to-ALP conversion probabil¬ 
ity is determined largely by the structure of the cluster 
magnetic field. 

Generically, clusters support magnetic fields of 0{fj,G) 
strength that are coherent over 1-10 kpc scales (for re¬ 
views see [21, 22]). For the Goma cluster, a detailed 
multi-scale model of the tangled magnetic field was con¬ 
structed in [23], and shown to give a good fit to observed 
Faraday rotation measures. 

For the Goma magnetic field of [23], we have simu¬ 
lated photon-to-ALP conversion of the ICM by numer¬ 
ically solving (2) along a set of sightlines that sample 
a given field-of-view. Along any given sightline, photons 
originate along each point in the cluster, and the net con¬ 
version probability is an average weighted by the thermal 
photon density at each point, which scales as ~ n^. 

The distorted X-ray spectrum is then given by. 


The intracluster medium (ICM) permeating galaxy 
clusters is a hot thermal plasma with temperatures, 
depending on the cluster, of between 2 and 10 keV. 
The cluster is visible in X-rays through the thermal 
bremsstrahlung of the ICM, with both continuum and 
line emission. 

The continuum emissivity is given by. 
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where A is a constant and the Gaunt factor g{i',T) is 
a slowly varying function of frequency. The strength of 
the atomic lines is set by both the overall metallicity 
and individual abundance of heavy elements, as well as 
the local ICM temperature, and can be calculated using 
programs such as AtomDB [17]. 

Thermal emission from the ICM has been extensively 
observed by several generations of X-ray satellites across 
the last few decades, providing an excellent fit to observa¬ 
tions. Figure 1 shows fits to emission from the Coma clus¬ 
ter (adapted from [18] and [19]), which hosts an approx¬ 
imately isothermal ICM with a temperature of 8.1 keV. 
Note the excellent quality of the fits, with residuals all 
below the 10% level. As a further illustration of the at¬ 
tained precision, we note that the possible 3.55 keV dark 


/distorted(^) — (1 (7/'—>.a) (^)) /intrinsic(^) ; (79) 

where {Pj^a){^) denotes the net conversion probability 
for the sightlines, averaged over the relevant field-of-view. 

For M < 10^^ GeV, the conversion probability satu¬ 
rates for the entire X-ray range, leading to a uniform 
reduction of the cluster luminosity to 2/3 of its orig¬ 
inal amount. However, as shown in equation 8, for 
precisely the most observationally interesting range of 
10^^ < M/GeV < 10 ^ 2 ^ the photon-to-ALP conversion 
probability is large at higher X-ray energies, while be¬ 
ing suppressed at lower energies. Thus, such photon-to- 
ALP conversion induces spectral distortions of the ther¬ 
mal ICM spectrum [24]. In Figure 2a we illustrate this by 
plotting {P^^a)i^) obtained from the simulations with 
nia = 0, M = 4 X 10^^ GeV and a field-of-view of size 
(100kpc)2 at the centre of the Coma cluster. 

Given the quality of the fits in Figure 1, the simulated 
photon-to-ALP conversion probabilities can be used to 
constrain the ALP parameters M and ma- By conserva¬ 
tively stipulating the absence of deviations in the spec¬ 
trum of 10% or greater, we find the constraints plotted in 
Figure 2b. For nia ^ 5 x 10 “^2 eV this provides a bound 
on the ALP coupling M <7 x 10^^ GeV, while for rria ^ 
5 X 10“^2 eV, the conversion probability is suppressed by 
the large ALP mass, and for uia ^ 1 x 10“^^ eV the cor- 
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(a) XMM-Newton best fit flux model from Coma, within (b) NuSTAR spectrum and best-fit single temperature model 
radius of ~ 300 kpc [18]. for the central 12' X 12' of the Coma cluster [19]. 

FIG. 1: Thermal X-ray spectrum from Coma as seen by XMM-Newton and NuSTAR. 


responding bound on the coupling becomes weaker than 
existing astrophysical constraints. 

We therefore conclude that constraints on spectral de¬ 
viations of the cluster thermal bremsstrahlung spectrum 
can provide the strongest astrophysical bounds to date on 
the ALP-photon coupling, excluding an interesting range 
of parameters that will be probed in upcoming laboratory 
experiments. 


Further properties of the distorted spectrum 

There are several further properties, and thus poten¬ 
tially correlated signals, of an ALP-distorted cluster spec¬ 
trum. 

First, along a single line of sight photon-to-ALP con¬ 
version is highly stochastic and energy-dependent. These 
features can provide several additional observational sig¬ 
natures, which are observable for small field-of-views. In 
Figure 3a we illustrate this by showing the averaged con¬ 
version probability for a (5kpc)^ field-of-view. Since the 
simulated field is random in nature, the exact energy de¬ 
pendence of the real field will differ, but the presence of 
these sinusoidal oscillations at X-ray energies is a charac¬ 
teristic feature of ALP-photon conversion, with the fre¬ 
quency of these oscillations decreasing with increasing 
energy. 

When averaging over a large field-of-view of size far 
greater than the coherence length of the cluster magnetic 
field (as in Figure 2a), these variations wash out. How¬ 
ever, by extracting the cluster X-ray spectrum across a 
small region it may be possible to directly observe such 
oscillations. The best candidates for this are very nearby 
clusters, where the fixed telescope angular resolution cor¬ 
responds to the smallest physical scales. 

For example, the Virgo cluster is at a redshift of z = 


0.004 and there is approximately 500 ks of observation 
time on Virgo with the Chandra telescope, whose superb 
imaging capabilities gives it arcsecond resolution. By 
considering spectra extracted from small angular regions 
instead of the full field-of-view, one could search for the 
presence (or absence) of such oscillatory features. 

There is another approach to searching for these small 
scale features. Clusters are large Mpc-scale objects, and 
gradients in their internal structure will generally also be 
large-scale. In the absence of photon-axion conversion, 
the bremsstrahlung photon count from nearby pixels on 
a detector is expected to be identical, subject to the vari¬ 
ations due to Poisson statistics. Photon-ALP conversion 
can lead to significant variations in conversion probabil¬ 
ities from nearby sightlines, and so in the presence of 
ALPs the photon count from nearby pixels should have 
a greater variation than would be expected simply from 
Poisson statistics. 

Furthermore, along a single line of sight the X-ray 
emission becomes highly polarised. In Figure 3b we 
show that with a small field-of-view this induced polari¬ 
sation would be observable, but again for larger field-of- 
views this effect gets washed out. The result is that for 
10^^ GeV ^ M < 10^^ GeV, the diffuse thermal ICM 
emission extracted from a small region would be po¬ 
larised. If an X-ray telescope with polarisation capability, 
such as the European LOFT or Ghinese XTP missions, is 
(finally, after forty years) flown, the presence or absence 
of polarised emission from the IGM could lead either to 
the detection of ALPs or strong bounds on their coupling. 
While optical polarisation induced by ALP-photon oscil¬ 
lations has been considered before (see for example [25]), 
the advantage of a galaxy cluster environment is efficient 
X-ray conversion and better studied magnetic fields. 

The small-scale spectral oscillations can also be sought 
by considering X-ray point sources that are fortuitously 
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(a) Averaged conversion probabilities for 
M = 4 X 10^^ GeV, = OeV, for a 100 kpc 
field-of-view. The dashed line represents the saturation 
value of 1/3. 
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(b) Astrophysical bounds from the lack of an observation 
of a gamma-ray burst with SN1987a in red, bounds on 
the parameter space from this study in blue. Because of 
uncertainties in the cluster magnetic field, the precise 
values of M excluded may be uncertain by a factor of 
two. 


FIG. 2: Averaged conversion photon-to-ALP conversion probability and constraints on M. 



(a) The averaged conversion probabilities. 


FIG. 3: Conversion probabilities and polarisation for 

and m; 



(b) The difference between conversion probabilities for the 
two polarisations of X-rays, showing the induced polarisation 
caused by ALP-photon oscillations. 

small 5 kpc field-of-view, with coupling M = 4 x 10^^ GeV 
1 TOq = 0 eV. 


observed through a cluster. Such a search has been car¬ 
ried out for the Hydra A cluster in [26] using a central 
AGN source. However this is not an ideal source. Hydra 
A is at a redshift oi z = 0.06, more than ten times further 
away than e.g. the Virgo cluster, and furthermore the 
AGN is optically thick, with significant intrinsic absorp¬ 
tion at lower X-ray energies, making it harder to separate 
photon-ALP leakage from this effect. 

Another interesting effect would be that the spectral 
distortions would have a radial dependence within a clus¬ 
ter. The free electron density decreases with radius, and 
observations indicate the magnetic held strength also de¬ 
creases with radial distance from the centre of the cluster. 
Photons travelling along off-centred sightlines will then 
convert at a reduced rate at higher energies (see equation 
6), and the ‘step’ in the conversion probability appear at 
lower energies, as implied by equation 8. This behaviour 


would be characteristic of ALP-induced distortions. 


Conclusions and Outlook 

The main point of this paper is that thermal emis¬ 
sion from galaxy clusters provides a well characterised 
X-ray source, arising within a magnetic held environ¬ 
ment that would lead to signihcant spectral distortion 
for ALP masses ma ^ I x 10“^^ eV and ALP-photon 
couplings 10^^ GeV ^ M < 7 x 10^^ GeV. As such spec¬ 
tral distortions do not appear to be observed, this places 
highly competitive, and potentially the most stringent, 
bounds on ALP parameter space. We should of course 
also caution that bounds on M can only ever be as good 
as knowledge of the astrophysical magnetic held, which 
for clusters is probably currently uncertain up to a factor 
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of two. 

We have used the absence of any reported large devi¬ 
ations from a thermal spectrum to bound the ALP pa¬ 
rameter space. Existing analyses however are not opti¬ 
mised for ALP searches. For these, it is beneficial to ex¬ 
tract spectra over a reduced field-of-view and search for 
energy-dependent oscillations. We plan to pursue this 
approach in future. 

An attractive reason to study X-ray emission from clus¬ 
ters and its connection to ALPs is that the data is only 
going to get better. The three current satellites XMM- 
Newton, Chandra, and Suzaku, all use CCD detectors 
with an intrinsic resolution of AA ~ 100 eV. Next year 
ASTRO-H will launch, which will allow a far greater pre¬ 
cision to be obtained due to the 7 eV resolution available 
from the microcalorimeters on board its Soft X-ray Spec¬ 
trometer. This improved resolution is highly relevant for 
searches for the sinusoidal spectral oscillations character¬ 
istic of ALPs. Over the longer term, ATHENA has been 
approved by ESA for a 2028 launch. 

Knowledge of cluster magnetic fields will also 
be enhanced over the short, medium and longer 
term. Through observing (or not) upscattered Inverse 
Compton-CMB photons, hard X-ray imaging telescopes 
will either produce direct measurements of the mag¬ 
netic field, or more stringent lower bounds on the field 
strength. Over the next decade, the Square Kilometer 
Array will also come online. With a large increase in the 
number of available radio sources, this will improve Fara¬ 
day rotation measurements of galaxy clusters’ magnetic 
fields. 
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